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We study design challenges associated with realizing a ground state quantum computer. In such 
a computer, the energy gap between the ground state and first excited state must be sufficiently 
large to prevent disruptive excitations. Here, an estimate is provided of this gap as a function of 
computer size. We then address the problem of detecting the output of a ground state quantum 
computer. It is shown that the exponential detection difficulties that appear to be present at first 
can be overcome in a straightforward manner by small design changes. 



I. INTRODUCTION 



Recently, there has been intense interest among researchers in the possibility of designing quantum computers 
that calculate using the remarkable properties of quantum mechanics |2|-|^] . Although the potential power of quantum 
computation algorithms is enticing, there is great difficulty associated with actually fabricating a quantum computer 
in the laboratory. A variety of schemes have been suggested pj-p^ and progress has been encouraging, but the 
feasibility of realizing a useful quantum computer is still unclear. 

In a recent paper, we proposed a novel "ground state quantum computation" approach that could circumvent some 
of the main problems with traditional quantum computer designs p9| . This approach replaces the progress of a 
usual, time-dependent quantum computation with a single, time-independent state. To see how this works, let us 
suppose that a quantum computation algorithm consists of the evolution of M qubits through N steps defined by 
2M ^ unitary matrices Uj, j = 1, . . . , N. To define the state of M qubits at one step of the algorithm requires 2*^ 
amplitudes. To describe the M qubits at every step of the algorithm, from before it begins until after it ends, requires 
(A^ + 1) X 2*^ amplitudes. If we do not demand that the qubits evolve simultaneously from step to step, allowing qubit 
#1 be at step 2 while qubit #3 is at step 6, then (2(A^-|- 1))^^ amplitudes are required to map out the development of 
the qubits. Let us suppose that we collect these (2(A^ -I- 1))*^ amplitudes into a state |^) defined on a Hilbert space 
of dimension {2{N + 1))*^. This state will contain all of the information in a time-dependent quantum computation, 
but the state itself will be completely time-independent. Instead of developing through time in accordance with an 
algorithm, the state will develop through Hilbert space in accordance with the algorithm. 

How can we explicitly describe this development through Hilbert space? The projection of |^) onto some 2*^ 
dimensional subspace will contain the 2*^ amplitudes necessary to describe the state of the M qubits when they are 
all at step and have undergone no unitary evolution. Let us call this projection Pq |^)- More generally, let us 
call Pj l^f) the projection onto the 2^ dimensional subspace that describes the state of the M qubits when they are 
at step J, J = 0, ... , N. Suppose that we define an operator Aj^Q that carries the 2^^ basis vectors of the subspace 
associated with Pq into the 2*^ basis vectors of the subspace associated with Pj. Then, |^) develops in accordance 
with a quantum computation algorithm provided that 

I*) = UjUj^i . . .UiA.^oPo I*) (1) 

for j = 1,...,N. 

This formal notion makes it possible to propose a new ground state approach to quantum computation. In ground 
state quantum computation, we do not make a register of qubits evolve in time by subjecting it to a series of time- 
dependent Hamiltonians. Instead, we perform calculations by manufacturing a Hamiltonian H whose ground state 
develops according to the equation (|l|). This is described precisely in reference In that paper, an appropriate 
Hamiltonian H is found. It is comprised of a sum of (i) one-body terms denoted like h'l{Ua,k) where Ua,k indicates the 
unitary evolution of qubit a at algorithmic step k of the calculation and (ii) two-body terms designated /i^ ^(CNOT) 
associated with a controlled not of qubit b by qubit a at step j. In [|9l, a possible physical realization of H is suggested 
using states localized on quantum dots to comprise the (2(7V -t- 1))^ dimensional Hilbert space. 

This ground state, time- independent approach has the attractive characteristic that it avoids traditional decoherence 
problems associated with time evolution. However, it does have its own challenging aspects that need to be addressed. 
Many of these challenges concern the scaling of a ground state quantum computer - the feasibility of making such a 
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computer larger and larger. In this paper, we address two of the most important considerations involved in increasing 
the size of a ground state quantum computer. First, we study how the energy gap between the ground state and 
first excited state of the Hamiltonian H will depend on N and M. Clearly, if a ground state quantum computer is to 
function properly the gap must be large enough in energy that the computer will reliably remain in its ground state. 
Second, wc investigate the problem of measuring the outcome of a ground state quantum computation. In its most 
naive form, a ground state quantum computer would become very difficult to probe as it grew in size. We propose 
several means of avoiding this difficulty. 



II. GAP 



If a ground state quantum computer is raised into an excited state, its wavefunction can no longer be relied upon to 
satisfy equation (|l|), and the computer therefore does not compute correctly. To avoid such excitations, the computer 
must possess a sufficiently large energy gap between ground state and first excited state. The gap must be significantly 
larger than the available thermal energy ksT, for example. 

As a ground state quantum computer grows in size, its gap will decrease, limiting the length of practical compu- 
tations. Here, we study this limit, describing the size dependence of the gap of the Hamiltonian described in [ p^ . 
We argue that the gap shrinks approximately like 1/{N +1)^ and prove the existence of a lower bound that scales as 
1/(^+1)^ 



A. Single Qubit 

To obtain these quantitative estimates of the gap, we first consider the case of a single qubit computer. Here, the 
Hamiltonian is simply H = '^f^i fi^{Ui) where 



h'(U) = e 



cj_iC\-i + c}a - {cjua^i + h.c.)] . (2) 



is associated with the development of the single qubit from step i — 1 to step i of the calculation. It is convenient 
to make a unitary transformation from the Ci operators to new operators (HJ^]^J7j)Ci. This changes the form of the 

Hamiltonian to H — X^il^i ^'(-^) where every Ui that appears in H has been replaced by the 2 by 2 identity matrix /. 
To determine the eigenspectrum of this new H , we solve the determinantal equation det(H — E) — 0. The determinant 
is evaluated by deriving and solving a recursion relation on matrices of increasing size. We find that 

detiH - E) ^ _ (3) 



where k = (1 ^ ^) + ^(1 — — 1. Setting the determinant (||) to zero yields eigenenergies E = Em = 2e(l — 
cos{-Km/2{N + 1))), where m is an integer between and 2N + 1. The ground state has energy E = Eq = 0, and 
the first excited state has energy E = Ei = 2e(l - cos{tt/2{N + 1))) ^ en^/{2{N + l))^ for large N. Thus, the gap 
decreases as l/(-/V +1)^. This is true for a single qubit, and also for any number of non- interacting qubits. 



B. One CNOT Gate 



Of course, a useful quantum computer must have interactions among qubits, so the behavior of the gap must be 
examined when interactions are present. To begin, we address the case of exactly two qubits interacting through 
exactly one CNOT gate. The full Hamiltonian includes single qubit h'^{Ui) terms and one W{CNOT) interaction 
term at stage j. 

We begin by examining the Hamiltonian with the (CNOT) term omitted. Without the (CNOT) term, the 
computer has two disjoint regions for each qubit, one "upstream" of the omitted CNOT gate consisting of stages 
to j — 1 and one "downstream" of the omitted CNOT gate consisting of stages j to N. Since an electron in one of 
the disjoint regions will possess the eigenspectrum of a single, non-interacting qubit, the first excited state in such a 
region will have an amount of energy of order l/(iV -I- 1)^. If we neglect such "high-energy" states, only the (doubly 
degenerate) ground states of the two regions make important contributions to the electronic state of each qubit. This 
means that each qubit has four available states, leading to an effective 16 dimensional Hilbert space for the two qubit 
system. 
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It is straightforward to diagonalize the interaction Hamiltonian h^{CNOT) analytically in this 16 x 16 basis. The 
result is a (fourfold degenerate) ground state of the computer with zero energy, an (eightfold degenerate) first excited 
state with energy e/(j){N — j + 1), and a (fourfold degenerate) second excited state with energy e/{N — j + 1)^ + 
+ So, the energy of the gap in the 16 dimensional Hilbert space scales as e/(j)(iV—j + l) ^ 1/7V^. 

What relationship does the gap in this 16 dimensional Hilbcrt space have to the exact gap of the system? The 
(fourfold degenerate) ground state in this 16 dimensional Hilbert space is, in fact, the exact ground state in the whole 
Hilbert space. Hence, the (eightfold degenerate) first excited state in the 16 dimensional Hilbert space is actually 
orthogonal to the exact ground state of the system. It follows that the quantity e/ {j){N — j + 1) represents a rigorous 
variational upper bound to the exact first excited state energy of the system. 

The variational upper bound should provide a reasonable estimate of the true value of the gap. However, for our 
purposes we are perhaps more interested in having a guaranteed lower bound to the gap. Such a lower bound would 
ensure that, when less than a specified amount of energy is available, the computer will not experience a disruptive 
excitation. As it turns out, it is possible to show that the gap has a rigorous lower bound of a/{N + 1)^ for some real 
positive a. The following is an argument by contradiction. 

We know the ground state of the Hamiltonian H has energy zero. Suppose that is some state of the 2-particle 
system that is orthogonal the ground state of H. Assume that the expectation value H j-i/)) satisfies 



\H\ 



<«/(iV+l)4^i?l, 



(4) 



To evaluate the left hand side and draw a contradiction, we split the Hamiltonian H into Hq, which consists of 
only single body terms, and Hi = (CNOT), which consists of only interaction terms. Both Hq and Hi, it is 
straightforward to show, are positive semi-definite. Consider the form of j?/;) in a basis of eigenstates of Hq 



n,i 



(5) 



where the i labels the degenerate eigenstates with the nth eigenenergy. Saying that \t/j) is orthogonal to the (fourfold 
degenerate) ground state of H essentially means that \tp) has no contribution from the four states \(t>n,i) for which 
Hq + Hi \<j>n.i) — 0. Therefore, \tp) can consist of a superposition of eigenstates of Hq with eigenenergies greater than 
zero and the 12 zero energy eigenstates of Hq that are orthogonal to the (fourfold degenerate) ground state of the 
system. 

Further reflection shows that [ip) cannot involve exclusively the 12 zero energy eigenstates of Hq because then 
{tp\H\ilj) would go like 1/(A^ + 1)^, as we showed above, violating the assumption (Q). The state must therefore 
possess some contributions from excited eigenstates of Hq. These states have eigenenergies of at least e7r^/(2(iV-|- 1))^, 
as we saw in our single qubit analysis, so assumption (^ Hmits the contribution from such states to J2n>o % < 
i?jQ^gj./(e7r^/(2(iV + 1))^). This limit exists even though Hi is present since Hi is positive semi-definite and cannot 
decrease the expectation value produced by Hq. Hence we find that 



{il\H\lJj) = C*^0,iC«=0,j {(t>n=aA\ Hq + Hi |(/),t^Oj) + C* jCmj ((/)„,i| + ^1 |0mj 

i.j n>0..'m>0.i.j 



> 



- (C,,:Cn=0,i (0m,j| -ffl \<l)n=0,j) + C*=o,»Cmj (0n=O,j| Hi jt/'mj)) 

m>0,'i,j 



n>0,i 



(2(iV + l)) 



n>0.i 



E 



Cfi—OA 



{N + iy 



V4(iV+ 1)2 - 16. 



e^V(2(A^+l))2 



E' 

m>0,'i 



(iV+l)2 



(N + iy 



(6) 
(7) 

(8) 
(9) 



here — /i/(A^ + 1)^ is the most negative value of 
(0), however, since 



i>m>o,i \ Hi \(j)n=oj). This last inequality contradicts the assumption 



[N+iy 




(10) 



A^) ^lower 
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provided that a is chosen to be sufBciently smalL This contradiction shows that the assumption (jj) is not vahd. Since 
the ground state has energy zero, and any state orthogonal to the ground state has energy at least ii'iower' have a 
lower bound -Eiower on the value of the gap. 



C. Arbitrary Computer 

It is straightforward to apply these results to the case of M qubits interacting via an arbitrary number of CNOT 
gates. First of all, the variational upper bound on the gap of e/{i){N — j + 1) still holds. This is because different 
CNOT Hamiltonians h-'^ ^^{CNOT) commute with one another, so we can treat each separately when diagonalizing in 
a basis of zero energy, M-particle eigenstates of iJo- 

The lower bound on the gap of order 1/(A^ + 1)"' also still holds, which we demonstrate in the following way. The 
Hamiltonian consists of Hq, that governs the single qubit development between CNOT gates and the CNOT gate 
terms j^{CNOT) themselves. We begin by dividing Hq into parts labelled {HoY^ ^ where the index j, a, h suggests 

proximity to the CNOT gate controlled by Hamiltonian ^^{CNOT). Let (-ffo)a consist of terms that control the 
single qubit development of qubit a between CNOT gate j, a, h and the previous CNOT experienced by qubit a, terms 
that control the single qubit development of qubit a between CNOT gate a, h and the next CNOT experienced by 
qubit a, terms that control the single qubit development of qubit b between CNOT gate j, a, h and the previous CNOT 
experienced by qubit h, and terms that control the single qubit development of qubit h between CNOT gate j, a, h 
and the next CNOT experienced by qubit b. By this definition of the (-ffo)^ b, we have Hq = \ ^ ft(-ffo)a h+ extra 
positive semi-definite one body terms that are near the first or last stages of^ the computer. 

Now, with this division described, we are in a position to demonstrate the lower bound of order l/(iV+l)"'. Suppose 
that some M-particle state jV') is orthogonal to the (2*^ fold degenerate) ground state of the system. It is possible to 
write \^) in the form 



IV') = '^Cn,i \(l>n,i/ 



(11) 



where the \(l)n,i) are Af -particle eigenstates of Hq. Because \4>) is orthogonal to the ground state of the system, each 
term \(j)n,i} that appears in lip) must satisfy (i?o)a & + fc(^^*^^) I0ri,i) 7^ for some CNOT gate j, a, b. Let us call 
V'a bj 1^*5 °f components c„_j; \(f>n,i) that satisfy (i?o)a bi^^^'^) \4'nA) 7^ 0. If any \(f)n,i) could belong 

ipi b), we include it in every possible ipi b/- Then, 



in more than one 



j,a,b 



j-a.b 



(.Ho)U + KbiCNOT) 



> 



{N + 1)4 - 2{N + 1)4 {N + 1)4 ' 



(12) 
(13) 



where we have made use of the lower bound a/ {N + 1)^ derived in the last section. The last inequality holds because 



every component of \ip) appears in at least one of the 



holds for an arbitrary number of qubits and CNOT gates. 



This result shows that a lower bound ~ i/{N + 1)^ 



III. DETECTION 



The task of detecting the result of a ground state computation seems daunting at first. Each qubit in the computer 
must be measured in the final stage, which it only visits with probability 1/{N + 1). Since there are M qubits, the 
probability of making a successful measurement scales as l/(iV-|- 1)^^. This problem is not, however, unsurmountable. 
In the following we describe several schemes that circumvent the detection scaling problems. 

In the case of Grover's algorithm, or any other algorithm for which the final state of the qubits factors, it is actually 
possible to guarantee a successful measurement. This is accomplished by placing an additional readout electron just 
after the final stage of each qubit, as in Fig. |l|. Since the qubit electron will be localized in one of the two dots 
in the final stage, the readout electron will tend to get pushed toward the other dot by the Coulomb interaction. 
In the ground state of the entire system, qubit electrons and readout electrons included, a readout electron will be 
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localized to the left (right) if its qubit is in the right (left) at the final stage. The outcome of the algorithm can 
always be determined by detecting the position of the readout electrons, even though the qubit electrons have only 
a small probability of residing at the final stage. In effect, the readout electrons are providing a continuous, passive 
measurement of the final state of the computer. Although this readout electron method is guaranteed to work for an 
algorithm for which the final state of the qubits factors, it will not work at all if the final state of the qubits does not 
factor. For a general algorithm, another approach is necessary. 

An approach which applies to arbitrary algorithms involves the adjusting the Hamiltonian at the final stage iV 
for each qubit a. Suppose that the operator Cl ^ is replaced by fiCl everywhere it appears in the Hamiltonian 
and Ca,N by pCa.N, where /? is a small fraction. All algorithms will still work just as before, but we are "tipping" 
the computer toward the final stage so that the qubits reside there more often. Then, it follows that in the ground 
state of the system, each qubit has times greater amplitude on the final stage than on the previous stages. The 
probability of detecting all qubits on the final stage is of order 1/(1 + /?^A^)*^. If (3 is set to be of order, say 1/vM]V, 
we find that the probability of all qubits being at the final stage goes as approximately (1 — , which approaches 

exp(— 1). It only takes two or three attempts to catch all of the qubits at the final stage. 

Of course, the change in the Hamiltonian will effect the gap. If the final operators are scaled by a factor /3 then the 
quantity det(-ff — E) of non-interacting qubits will change to 

det(i/ - E) = ,^(-+i)^i_li^(fc2(iv+r) _ ^ _ i)(,...r _ 

Setting this determinant to zero, we find that the gap of non-interacting qubits still scales roughly as 1/(A'' + 1)^ for 
arbitrary [3 between zero and one. Once CNOT gates arc included, however, the gap will have a variational upper 
bound of order 1/{N + l)(iV + l//3^) and a lower bound of Ei^^^j. = a/{N + lf{N + l/(3'^ f. If /? = 1/VMN, then 
the upper bound is 1/{NM + N){N + 1) and the lower bound E'lower = "/(-^ + 1)^(^M + N)^. 

Another technique for alleviating measurement problems is to "synchronize" the arrival of the qubits at the final 
stage. In our CNOT gate, the target electron cannot proceed beyond the gate until the control electron has. A 
controlled "identity" (CID) gate could be constructed that would function similarly, preventing a target electron from 
proceeding beyond the gate until after a control electron, but always subjecting the target qubit to an identity operation 
and never a NOT operation. With this gate, the arrival of the qubits at the final stage could be "synchronized." 
Suppose that each qubit controls the entry of the next qubit to a ground state quantum computer's final stages, 
using a CID gate. Then, whenever qubit M is found in the final stage, all qubits are there. This could be useful for 
detection schemes, although it would not enhance the overall probability of finding the qubits at the final stage of the 
computer. 

IV. CONCLUSION 

In this paper, we have explored some important challenges to constructing a ground state quantum computer. 

First, we found upper and lower bounds for the energy gap between the computer's ground state and first excited 
state. The bounds provide guidelines to making a computer of a specified size that can be relied upon to remain in 
its ground state. Next, several schemes were presented for easing qubit detection difficulties. These schemes indicate 
how to probe a ground state quantum computer so that it will yield output with certainty or at least high probability. 
It is hoped that the analysis of energy gap and detection in this paper complements our initial proposal and eases the 
task of designing and fabricating a ground state quantum computer in the laboratory. 
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FIG. 1. Additional readout electrons are placed after the final stage of each qubit. The position of these additional electrons 
indicates the outcome of the algorithm. 
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